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DYNAMIC STABILITY OF A 

4.6-METER-DIAMETER 120° CONICAL SPACECRAFT AT 
MACH NUMBERS FROM 0.78 TO 0.48 IN A SIMULATED 
MARTIAN ENVIRONMENT 

By Charles H. Whitlock and Richard J. Bendura 
Langley Research Center 

SUMMARY 

Subsonic dynamic stability characteristics of a large-size 120° conical spacecraft 
have been determined from flight data in a simulated Martian environment. The results 
cover a Mach number range from 0.78 to 0.48, and a Reynolds number range from 
200 000 to 100 000 (based on 15-foot (4.6 meter) spacecraft diameter). Assumptions of 
the analysis include linear static aerodynamics and invariant dynamic stability derivative 
Cm Q + Cm^, with angle of attack. Reduction of the data on a cycle -by-cycle basis gives 

Cm a + Cm^ values near -0.135 with an uncertainty of ±0.105 over the range of the test. 

Trajectory simulation studies indicate that a reasonable simulation can be obtained by 
using a damping coefficient of -0.135, and that ±0.105 is a realistic value for the accuracy 
limit. 


INTRODUCTION 

Large-angle blunted cones are being considered as configurations for entry into 
planetary atmospheres. A blunt 120° total angle conical shape is proposed for several 
future programs and was used for the Planetary Entry Parachute Program. The dynamic 
stability characteristics of this shape are under investigation in ground facilities. Ref- 
erences 1 and 2 present dynamic stability derivative (C m q + C m ^ data for both low sub- 
sonic and hypersonic velocities on small-scale models. The purpose of this report is to 
present the dynamic stability characteristics of a large -size 120° conical spacecraft 
determined from free -flight data at Mach number and Reynolds number conditions that 
simulate a possible Martian environment. 

The results were obtained from the first flight test of the balloon-launched series 
of the Planetary Entry Parachute Program. Details of the flight are described in refer- 
ences 3 and 4. By utilizing the spacecraft attitude histories obtained in reference 5, 



values for the dynamic stability derivative C m q + C m( j are obtained at Mach numbers 
from 0.78 to 0.48 at earth altitudes between 138 000 feet (42.1 kilometers) and 143 000 feet 
(43.6 kilometers). Reynolds number (based on the 15-foot (4.6 meter) spacecraft diam- 
eter) varied from 200 000 to 100 000. 
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transformation matrix defined in equation (1) 


axial-force coefficient (positive rearward), 
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pitching-moment coefficient, 
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pitching -moment -curve slope, 
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da 


, per radian 
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, per radian 


normal-force coefficient, 


Normal force 


q’S 


9C N 


normal -force -curve slope, — — , per radian 

aa 

maximum body diameter, feet (meters) 

moment of inertia about pitch or yaw axis, slugs -feet^ (kilograms -meters^) 
moment of inertia about roll axis, slugs-feet^ (kilograms-meters2) 
mass of spacecraft, slugs (kilograms) 

roll, pitch, and yaw velocities, respectively, radians/second 

free -stream dynamic pressure, pounds force/foot^ (newtons/meter^) 

maximum cross-sectional area of body, feet 2 (meters^) 


k 


2 


t 


time from spacecraft separation, seconds 


I 

t' time to one-half amplitude of corrected angle-of-attack ratio squared, seconds 

u,v,w component of spacecraft linear velocity relative to earth along X-, Y-, and 

Z-axis, respectively, feet/second (meters/second) 

u'jV^w' component of spacecraft linear velocity relative to wind along X-, Y-, and 
Z-axis, respectively, feet/second (meters/second) 

V' free-stream velocity, feet/second (meters/second) 

V w xE’^w YE’^w ze w ind velocity component along Xg-, Yg-, and Zg-axis, 

respectively 

X,Y,Z body axes of spacecraft 

Xe,Ye,Ze earth-fixed axes 

Xg,YE>ZE spacecraft velocity components along Xg-, Yg-, and Zg-axis, respectively 
x C g center -of -gravity distance rearward from apex, feet (meters) 

a, 1 3 angles of attack in pitch and yaw, respectively, degrees 

rj resultant angle of attack, degrees 

9,<p,\l > pitch, roll, and yaw angles of spacecraft relative to earth-fixed axes, 

respectively, radians or degrees 

p atmospheric density, slugs/foot^ (kilograms/meter^) 

7 p flight-path angle in pitch, degrees 

Yy flight -path angle in yaw, degrees 

X Q damping parameter, 

Dots over symbols indicate derivatives with respect to time. 
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CONFIGURATION AND TEST ENVIRONMENT 

A sketch of the spacecraft configuration utilized for the first flight of the balloon- 
launched series of the Planetary Entry Parachute Program is shown in figure 1. The 
mission profile is shown in figure 2. The results described herein were obtained during 
the aeroshell coast portion of flight which occurred after parachute deployment and sub- 
sequent aeroshell separation from the payload. The center -of -gravity location shown in 
figure 1 is that after aeroshell separation. Mass properties are as follows: 

21.68 slugs (316.4 kg) 

x C g/D, rearward from apex 0.218 

lx 203 slugs-ft^ (275 kg-m^) 

I 159 slugs-ft^ (215 kg-m^) 

Mass and center of gravity are measured values. Roll inertia was estimated, and pitch - 
yaw inertia is an average value deduced from observed motion frequencies and wind- 
tunnel values for C mQ; . The Euler angle and aerodynamic angle system used in this 
analysis are those of reference 6. Figure 3 is a sketch showing the two systems. 

In terms of flight time, aeroshell separation occurred at approximately 14 seconds 
from spacecraft release and attitude data were obtained from onboard motion-picture 
cameras until 27.3 seconds of the release. Beyond this time, sight of the earth was lost 
by the rearward-facing cameras as the aeroshell approached apogee. Time histories of 
the Euler angles (0, 4/, and <p) are shown in figures 4, 5, and 6 (from ref. 5). Figure 4 

shows the 0 time history obtained from both cameras 1 and 2, and the ^ an d $ time 
histories are shown in figures 5 and 6, respectively. The Mach number and dynamic 
pressure histories during the data period are shown in figure 7. It should be noted that 
the dynamic pressure decreased (from 2.1 pounds per square foot (l00.5 N/m^) initially) 
by a factor of 3.28 during the data period. Reynolds number varied from 200 000 to 
100 000 based on spacecraft diameter. 

ANALYSIS AND RESULTS 

From the results of reference 5, values for the dynamic stability derivative 
Cniq + C m ^ are determined in the data-reduction section which follows. The validity 
of the results is tested by use of trajectory simulation techniques in a subsequent section. 

Data Reduction 

If the spacecraft attitude angles 0, <p, and ^ are known, free-stream velocity 
components along the body X-, Y-, and Z-axes can be calculated by use of the following 
equation: 
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(cos 6 sin 4d (-sin 9) 

(cos p cos <p + sin 4 / sin 9 sin tp) (cos 9 sin p) 

(sin p sin 9 cos p - cos p sin p) (cos 9 cos p) 


The earth-relative velocities X E , Y E , and Z E were determined from ground-track 
data. Wind velocities were obtained by means of an Arcasonde sounding rocket which 
was launched 1 hour after spacecraft release. Location of the Arcasonde during descent 
was within 20 miles of the aeroshell during the data period. It was assumed that the 
winds measured by the Arcasonde were those that existed during aeroshell flight. The 
wind velocity in the vertical direction V w gE was assumed to be zero. If the free- 
stream velocity components are known, spacecraft angles of attack can be calculated by 
use of the following relations (see fig. 3): 

a = tan~l ~ (2) 

u' 

(3 = tan -1 (3) 

u' 


rj = tan 


-i 


+ w 1 ' 


U' 


(4) 


Figure 8 shows the 77 time history, and figure 9 shows the a.,\ 3 crossplot for each 
camera. The fairing of curves between the data points in figure 8 was accomplished by 
utilizing the cross plots in figure 9. Average envelopes enclosing the 77 histories were 
also faired. The first peak on the camera 1 results was ignored because based on a 
detailed review of data used in reference 5, the point at 15.56 seconds was believed to 
contain greater than average error. Differences between the values obtained by the two 
cameras are attributed primarily to two sources of error. First, the basic method of 
obtaining the Euler angles from photographs is subject to error (see ref. 5), and second, 
it was established that the two cameras were misalined relative to each other and to the 
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longitudinal axis of the spacecraft. Camera misalinement would show up as a trim angle 
on the a, 1 3 crossplot. Examination of figure 9 shows camera 1 yawed relative to cam- 
era 2, but it is not known which camera is more closely alined with either the aeroshell 
trim or longitudinal axis (which may not coincide). However, based on figure 9, the mis- 
alinements between either camera and the trim axis appear to be small (less than 5°). 

In order to establish whether the spacecraft was stable or unstable, it was neces- 
sary to consider the general nature of the rj history. It has been generally observed 
(ref. 7, for example) that the amplitude ratio of the oscillatory motion varies inversely 
as the fourth root of the dynamic pressure ratio if C mQ! does not vary with Mach num- 
ber, the spacecraft roll rate is small, and there is no aerodynamic damping. Static sta- 
bility data derived from wind-tunnel results (shown in fig. 10) and the roll history in 
figure 6 indicate that the first two assumptions are reasonable in this case. Also a tra- 
jectory simulation of the test verified that the resultant angle -of -attack ratio did follow 
the inverse of the fourth root of the dynamic pressure ratio if no damping is present. As 
noted previously, both the dynamic pressure and Mach number varied during the data 
period. If the spacecraft contained no aerodynamic damping, the amplitude of rj would 
have increased by a factor of 1.33 during the time period 14.3 to 25.9 seconds based on 
the dynamic pressure history. The camera 1 results would have an amplitude of 32°, and 
the camera 2 data would show a peak of 36.6° at 25.9 seconds if no damping were present. 
Since neither amplitude increased by the factor of 1.33 (fourth root of the dynamic pres- 
sure ratio), it is concluded that the spacecraft did contain positive damping (negative 
Cmq + C m ^) during the data period. 

To evaluate the magnitude of C mQ + C m( j, the methods of reference 8 were used 
after a correction to account for the variable dynamic pressure was applied. Assumptions 
of the analysis include linear static aerodynamics and constant C ni q + C m ^ with angle of 

attack. Also the effects of trim were neglected. For each cycle the measured resultant 
angle -of -attack ratio squared was corrected by the square root of the inverse of the 
dynamic pressure ratio, and the damping parameter X Q was calculated by use of the 
relation 


A 


O 


0.693 

2t’ 


(5) 


where t' is the time to one -half amplitude of the corrected resultant angle -of -attack 
ratio squared (taken from similar log plots as shown in ref. 8). Damping derivatives 
were calculated for each cycle at an average Mach number from the relation 
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as well as those of 


The damping observed is assumed to include the effects of C m ^ 

C m q- Average values during each cycle for the parameter 4m/pSV' were used. The 
results of these calculations are shown in figure 11 at the average Mach number for each 
cycle. Also shown are the estimated uncertainty based on extreme fairings of the 
?7 envelopes, the uncertainty in inertia properties, and the possible errors in dynamic 
pressure values. Wind velocities were to the order of 120 feet per second (36.6 meters 
per second), and an uncertainty of 20 feet per second (6.1 meters per second) must be 
considered in the final reduced values. The estimated uncertainty in C m q + 
appears large when considered in terms of magnitude, but when compared with ballistics 
range results shown in reference 2, the flight test uncertainties appear reasonable. 

Trajectory Simulation 

In order to verify the general amplitude of the calculated C mq + C md , values, the 
flight trajectory was simulated on an IBM 7094 electronic data processing machine by 
using the six-degrees-of -freedom program described in reference 6. For brevity, the 
results from only camera 1 were considered. It was realized that a precise simulation 
would not be possible because of uncertainties in the spacecraft inertial properties and 
wind conditions used as program inputs. Also the Euler angle histories being simu- 
lated were subject to some error. (See ref. 5.) By considering the general level of 
Cm q + C m( j , however, it was believed that the general nature of the motion could be 
simulated. By use of outside bounds for + C m( j, it was expected that the estimated 

uncertainty values could also be substantiated . 

Review of the data in figure 11 indicates that values of -0.03 and -0.24 represent 
outside bounds of C m + C m . if Mach number effects are neglected. The average of 
these two values, -0.135, was selected for use in initial simulations to verify the general 
level of C m q + C m ^. Mass properties and static aerodynamics (fig. 10) discussed pre- 
viously were used. The initial conditions required were: (1) spacecraft position (Xg, 

Y E , and Z E ^; (2) spacecraft attitude (0, cf), and \p); (3) spacecraft velocity components 
relative to earth along body X-, Y-, and Z-axes (u, v, and w); and (4) spacecraft roll, 
pitch, and yaw velocities (p, q, and r). At an initial time of 14.17 seconds, spacecraft 
position was obtained from ground-tracking data, and attitude was taken from figures 4 
to 6. Velocity components along the body axes relative to Earth were calculated from 
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Initial roll, pitch, and yaw velocities were obtained from the relations: 

p = 0 - 0 sin 9 ^ 

q = 9 cos 0 + 0 cos 9 sin 0 \ (7) 

r = -9 sin 0 + 0 cos 9 cos 0 

The simulation process consisted of successive iteration to determine the combi- 
nation of initial conditions and aerodynamic trim angle which provided the best simulation 
of the flight Euler angle histories. After initial calculations with equation (7), the values 
for p, q, and r were adjusted slightly to provide a good simulation for the first- 
quarter cycle of motion (before damping could be effective). Then the effects of various 
trim angles were investigated until a combination which provided motions that were cen- 
tered about and in phase with the flight 9 and 0 histories was found. (Some readjust- 
ment of p, q, and r was necessary to simulate the first quarter cycle.) Finally, the 
initial value of C m q + C m ^ (-0.135) was changed first to -0.03 and then to -0.24 to 
determine the effects of the estimated accuracy limits on the spacecraft motions. 

The final initial conditions used in the simulation are: 


t, seconds 14.17 

Xe, feet (meters) 0 ( 0) 

Ye, feet (meters) 0 ( 0) 

Z E , feet (meters) -136150 (-41499) 

9 , radians 0.488 

0, radians 4.42 

0, radians 2.78 

u, feet per second (meters per second) 691 (211) 

v, feet per second (meters per second) 333 (101) 

w, feet per second (meters per second) -23 ( -7) 

p, radians per second 0.26 

q, radians per second 0.25 

r, radians per second -0.35 


The aerodynamic trim angles which gave the best overall results were 


^rim - 0-6° 
^trim ~ "1-9° 
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The final values for p, q, and r were changes from 0.312 rad/sec, 

0.912 rad/sec, and -0.298 rad/sec which were the respective values used in initial 
simulations. Changes were necessary because of inaccuracies caused by taking slopes 
of 9, <p, and for use in equations (7); and by the fact that each time the trim angle 

was varied, q and r had to be readjusted to simulate the first-quarter cycle of motion. 

The values of 0, p, and (p from the best simulation with C m q + C m( j = -0.135 
are compared with the flight data in figure 12. A precise simulation could not be obtained 
with constant aerodynamics and mass properties. The low-amplitude oscillations of 6 
and \p are not duplicated. The large-amplitude motions are in good agreement; hence, 
the 7] history is expected to be a reasonable simulation of the flight conditions. Com- 
parison of the simulation and flight values of 77 in figure 13 substantiates this conclusion. 
The inability to simulate the low-amplitude motions is probably caused by errors in the 
inertia values used, the possibility of nonsymmetric pitch and yaw inertias, and the proba- 
bility of products of inertia in the flight vehicle. Also the previously discussed uncer- 
tainty in the wind profile used in the simulation is a source of error. The a, ft crossplot 
from the simulation is shown in figure 14. The width of the loops and the low-amplitude 
values of 77 do not agree with the camera 1 results in figure 9. Comparison of the 
flight and simulation motion patterns with the typical motion patterns (below roll reson- 
ance) in reference 8 suggests that slight errors in the inertia properties are the most 
probable cause for the inability to simulate the motion patterns. Error in the calculated 
values of C m + C m • due to inertia uncertainties is expected to be small, however, 

4 

because the oscillation frequencies correlate and indicate that the pitch and yaw inertias 
are in the right range. If figure 13 is compared with figure 8 , it appears that the overall 
trend of the simulated 77 history is in good agreement with the motion characteristics 
deduced from both sets of camera data. It is concluded that a reasonable simulation of 
the flight motions could be obtained with C m q + C m ^ = -0.135. 

Simulated trajectories were also computed by using C m + C m • equal to -0.03 

4 

and -0.24 in order to test the accuracy limits established during data reduction. The 
same initial conditions and trim angles as described previously were used so that the 
first -quarter cycle of motion was simulated. The results of these computations are 
shown in figures 15 and 16. Figure 15 shows 8, ip, and (p histories for both values 
of C m q + C m ^ along with the flight data points. Variation of C m g + C m ^ does not 
significantly affect the low-amplitude values of 6 or ip which were not well simulated 
in the first place. In general, minimum damping appears to cause excessive amplitudes, 
and maximum Cm^ + C m ^ causes the simulation to begin to peak beneath the high- 
amplitude values during the latter portion of the data period. These same trends are 
evident from the comparison of the 77 histories shown in figure 16. The low-damping 
simulation causes excessive divergence, and the high-damping case maintains nearly 
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constant amplitude. The best fairings of flight data in figure 8 show a slight divergence; 
hence, it is concluded that -0.03 and -0.24 are reasonable accuracy limits. 


CONCLUDING REMARKS 

Dynamic stability derivatives were obtained from the flight test of a 4.6-meter- 
diameter 120° conical spacecraft in a simulated Martian environment. Comparison of 
the rate of divergence of the resultant angle-of-attack rj histories with predictions 
based on the dynamic pressure history indicate positive damping during the data period. 
Actual reduction of the data on a cycle -by- cycle basis by using faired rj envelopes 
gives dynamic stability derivative (C m q + C m( jJ values near -0.135 in the Mach number 
range from 0.78 to 0.48. An estimated uncertainty of ±0.105 is based on extreme fairings 
of the rj history, possible errors in the inertia values, and inaccuracies in the wind 
data. Trajectory simulation studies indicate that a reasonable simulation can be obtained 
by using a damping coefficient of -0.135 and that ±0.105 is a realistic value for the accu- 
racy limit. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 19, 1968, 

124-07-03-05-23. 
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Figure 1.- Side- and rear-view sketch of spacecraft. All dimensions and stations are in inches (cm). 








Figure 2.- Mission profile. 
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Figure 11.- Variation of C m „ + C mrV with Mach number 
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Figure 16.- Simulated n time histories using Cmq + C m ,j accuracy limits. 
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